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ABSTRACT 

We present constraints on the stellar initial mass function (IMF) in two ultra-faint dwarf (UFD) 
galaxies, Hercules and Leo IV, based on deep Hubble Space Telescope (HST) Advanced Camera for 
Surveys (ACS) imaging. The Hercules and Leo IV galaxies are extremely low luminosity [My — —6.2, 
—5.5), metal-poor (([Fe/H])= —2.4, —2.5) systems that have old stellar populations (> llGyr). 
Because they have long relaxation times, we can directly measure the low-mass stellar IMF by counting 
stars below the main-sequence turnofF without correcting for dynamical evolution. Over the stellar 
mass range probed by our data, 0.52 — 0.77 M0, the IMF is best fit by a power-law slope of a = 1.21q'5 
for Hercules and a = 1.3±0.8 for Leo IV. For Hercules, the IMF slope is more shallow than a Salpeter 
IMF (a = 2.35) at the 5.8-cr level, and a Kroupa IMF [a = 2.3 above 0.5 Afg) sat 5.4-cr level. We 

for Hercules, and OAltlf^ for 
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simultaneously fit for the binary fraction, finding /binary — 0-47_q ^4 
Leo IV. The UFD binary fractions are consistent with that inferred for Milky Way stars in the same 
mass range, despite very different metallicities. In contrast, the IMF slopes in the UFDs are shallower 
than other galactic environments. In the mass range 0.5 - 0.8 M©, we see a trend across the handful 
of galaxies with directly measured IMFs such that the power-law slopes become shallower (more 
bottom-light) with decreasing galactic velocity dispersion and metallicity. This trend is qualitatively 
consistent with results in elliptical galaxies inferred via indirect methods and is direct evidence for 
IMF variations with galactic environment. 

Subject headings: Local Group galaxies: dwarf galaxies: photometry galaxies: evolution galaxies: 
formation galaxies: stellar content 
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1. INTRODUCTION 

The stellar initial mass function (IMF) parameterizes 
the relative number of stars formed in a single age pop- 
ulation as a function of stellar mass. The IMF is fun- 
damental to all calculations of star formation rates and 
galaxy stellar masse s (for reviews seelBastian et al.ll2010l : 
iKroupaet aD[20TTI ). The classic [SalEeted (11955D IMF is 
a single power law with slope a = 2.35 {dN /dm (x to~", 
where N is the number of stars of mass m). Salpeter's 
original IMF was based on stars down to 0.4 Mq, and 
modern Milky Way studies i ndicate a break in the IMF 
slope below this mass scale. iKroupal ()2002( ) parameter- 
ized the Milky Way IMF as a broken power law, with 
a = 2.3 above 0.5 M© and a shallower a = 1.3 slope 

^ Based on observations made with the NASA/ESA Hubble 
Space Telescope, obtained at STScI, which is operated by AURA, 
Inc., under NASA contract NAS 5-26555. 

^ Astronomy Department, Yale University, New Haven, CT 
06520, USA; marla.geha@yale.edu 

^ Space Telescope Science Institute, 3700 San Martin 
Drive, Baltimore, MD 21218, USA; tbrown@stsci.edu, tumlin- 
son@stsci.edu 

* University of Physics &; Astronomy, University of California 
Irvine, 4129 Frederick Reines Hall, Irvine, CA 92697, USA 

^ Center for Galaxy Evolution Fellow 

® Department of Physics and Astronomy, University of Victo- 
ria, P.O. Box 3055, Victoria, BC, V8W 3P6, Canada 

"^ Departamento de Astronomia, Universidad de Chile, Casilla 
36-D, Santiago, Chile 

* Observatories of the Carnegie Institution of Washington, 813 
Santa Barbara Street, Pasadena, CA 91101, USA 

® UCO/Lick Observatory and Department of Astronomy and 
Astrophysics, University of California, Santa Cruz, CA 95064, 
USA 



below this mass. iChabrieil (|2003l |2005[ ) parameterized 
the IMF below 1 Mq as a log-normal distribution with a 
characteristic mass rric ~O.2M0. 

One may expect that the low mass IMF should depend 
on the physical properties of the stellar birth cloud such 
as the gas density, metallicity, or turbulent velocity (e.g.,; 
Larson 2005; Bate 2009'; iMvers et all 12011 '; "Marks_et_al] 
2012: Hopkins .2012.), however, IMF observations are 



large l y invariant within t he Milky Way (jBochanski et al.l 
I2OIOI : ICovev et afl 120081 ). Recent indirect studies sug- 
gest that the low mass IMF slope does vary outside 
the Milky Way and may be a function of the global 
galactic potential: studies of integrated line strengths, 
kinematics and gravitational lensing studies of elliptical 
galaxies appear to favor IMFs that become increasingly 
bottom-heavy (with IMF slopes similar to or steeper than 
Salpeter) toward higher galaxy masses ("T reu et"aLll2010 ; 
|Van Dokkuni & Conroy 2011; Cappcllar i et al.l 12012a : 
iDutton et al.ll2012l: IConr ov & van Dokkum||20il). 

The ultra-faint dwarf (UFD) galaxies are a recent class 
of diffuse Galactic satellites disco vered in Sloan Digi- 
tal Sky Survey (SDSS ) data ('e.g.. I WiUman et al.l [20051 : 
iBelokurov et al.ll2007l ) . The UFDs are the l east luminous 
(-8 < Mv < -1.5: iMartin et al.l 12008; 'Mufioz et al" 



201C ) and most dark matter-dominated (Simon fc Geha 



20071 ) galaxies known. Their a verage metallicities are less 



than typical globular clusters (jKirby et al.ll2008l ). Anal- 
ysis of HST photometry implies stellar ag es at least as 
old a s the oldest known globular clusters (|Brown et al.l 
l2012f ). Thus, the UFDs are a prime environment to 
test predicted IMF variations with the temperature, den- 



sity, or cosmic e poch of the star-forniing environment 
(jTiimlinsonl [20071 ) ■ 

Direct estimates of the IMF based on counting re- 
solved main-sequence stars are largely l imited to the 
nearb y Galactic field and star clusters (jBastian et alJ 
l2010f l. Milky Way open clusters provide a relatively 
large stellar mass range ov er which to measu re the 
IMF (0.08 - 7 Mq, e.g., iMoraux et all [200l . but 
are more metal-rich than the UFDs. At metallici- 
ties less than [Fe/H] < —1, HST studies of Galac- 
tic globular clusters probe the IMF down to main se - 
quence masses between 0.2 — 0.7 Mq (|Paust et al.ll2010l ). 
However, dynamical evolution, such as mass segre- 
gation and evaporation, can significantly change the 
slope of the m ass function ([Vesperini fc Heggid 119971 : 
iBaumgardt fc M akino 2003). Interestingly, the observed 
MF slop es correlate contrary to expectation with concen- 
tration (|De Marchi et al.]l2010( ). The observed correla- 
tion can be understood as dynamical evoluti on combined 
with either gas expulsion of residual gas ([Marks et al.l 
|2012| ). or related to orbital properties and the degree of 
tidal stripping (Kiipper et al. in prep.). 

Dwarf galaxies have metallicities similar to or lower 
than Galactic globular clusters, but have relaxation times 
longer than a Hubble time and therefore do not require 
corrections for dynamical evolution. Dwarf galaxies in 
which the low mass IMF has been directly measured 
are the Milky Way satellites Ursa Minor {My = -9.2, 
[Fe/H] = -2.0), Draco (My = -8.6, [Fe/H] = -2.0), 
and the SmaU Magellanic Cloud (SMC; My = -15; 
[Fe/H] = -1.2). IWvse et al.l (l200l used HST/WFPC2 
data to conclude that the Ursa Minor IMF is consistent 
with a power law slope a = 1.8 o ver the mass range 
0.4 - 0.7 Mq. iGrillmair et al.l (|1998f ) found a power law 
slope for the Draco dwarf galaxy between 2.1 < a < 2.3 
for an assumed age of 12Gyr based on ggr/W FPC2 
imaging extending to 0.6 Mq. iKalirai et al.l (|2013D used 
HST/ ACS data to conclude that the IMF of the SMC 
has a power law slope of a = 1.90lg'}5 '^^^^ ^^^ mass 
range 0.37- 0.93 Mq. 

Using the HST/ ACS we are undertaking a deep imag- 
ing survey of UFDs reaching several magnitudes below 
the main-sequence turnoff. The program includes the 
Milky Way satellites Hercules, Leo IV, Ursa Major I, 
Bootes I, C oma Berenices, and Canes Venatici II. In 
iBrown et al.l (|2012f ). we presented a preliminary analy- 
sis of the stellar populations in the first three galaxies 
for which data had been taken. We concluded that all 
stars in these galaxies were older than 11 Gyr and that 
star formation lasted less than 2 Gyr. 

Here we present a companion analysis of the IMFs for 
Leo IV and Hercules. Althou gh Ursa Maj o r I w as in- 
cluded in the age analysis of IBrown et al.l (|2012l ). the 
ACS catalog of this galaxy has significantly fewer stars, 
and these stars were observed over a much wider area 
(specifically, 9 ACS tiles in Ursa Major I, compared to 1 
and 2 tiles in Leo IV and Hercules, respectively). While 
Ursa Major I is closer than either Leo IV or Hercules, 
it is subject to significant contamination from image ar- 
tifacts and field contamination, and does not provide a 
good constraint on the IMF. We therefore do not include 
Ursa Major I in our analysis. 

In §[2]we describe the HST/ ACS data and construction 



TABLE 1 
UFD Galaxy Properties and IMF Results 



Row Quantity 



Units 



Hercules 



Leo IV 



(1) 


a (J2000) 


h : m : s 


16:31:05 


11:32:57 


(2) 


<5 (J2000) 


o . / . // 


+ 12:47:18 


-00:31:00 


(3) 


(m-M)v 


mag 


20.90 


21.15 


(4) 


Distance 


kpc 


135 


156 


(5) 


E(B-V) 


mag 


0.08 


0.06 


(6) 


Mv 


mag 


-6.2 ±0.4 


-5.5 ±0.3 


(7) 


Lv 


Lq 


2.6tl-l X 10" 


1.41°-* X 10* 


(8) 


'r-cB 


pc 


230 


130 


(9) 


a 


km s 1 


5.1 ±0.9 


3.3 ±1.7 


(10) 


([Fe/H]> 


dex 


-2.41 


-2.54 


(11) 


^rclax 


years 


3 X 1012 


2 X 1012 


IMF Results 








(12) 


Mass range 


Mq 


0.52 - 0.76 


0.54 - 0.77 


(13) 


A^star 




2350 


1054 


(14) 


a 




1 2+°-" 
1-^-0.5 


1.3 ±0.8 


(15) 


rric 


Mq 


n 4+0-9 


4+2-1 


(16) 


/binary 




r, ,0+0.20 

0-48_o.i2 


0.47«:?? 



take 



from IMarti] 



Not e. — (1-2) Right ascension and dcclinat 
IliOO^ )- (3-5) The distance modulus, distance and reddening are best fitti 
values from the H ST C MDs. (5-7) Absolu te magnitudes and effective radii from 
Hjand et aTl iM6ii ) and lijand et al.l |56i(l| ) for Hercules a nd Leo IV, respe c tively . 
(7-8 ) Velocity d i spersio n and average metallicity from lijlmon fc (Ijehal jliOOl) 
and IKirbv et al.l |2U11| ). (9) The two-b ody relaxation is ca l culate d using the 
quantities in (5-7) and Equation 1.38 in liJinnev ^ TremaTTT^ HUU^ . (10) Stel- 
lar mass range (10) and number of stars (11) included in our IMF analysis. 



(11-13) Estimated IMF parameter, 
fractions (14) are given for the po 



1 the 
la 



■ 0.52 - 0.77 M, 



■ fits 



0- 



of the IMFs. In §[3l we use these data to constrain the 
IMF slope and binary fraction over the stellar mass range 
0.52 — 0.77 Mq. In §|31 we discuss implications of our 
IMFs for the UFDs. In §[5l we compare these results to 
other direct IMF measurements and discuss them in a 
broader cosmological context. 

2. OBSERVATIONS AND DATA REDUCTION 

The Hercules and L eo IV UFDs were discovered by 
iBelokurov et al.l (|2007D as statistically significant over- 
densities of stars in the SDSS. As listed in Table 1, Her- 
cules and Leo IV have average metallicities of ([Fe/H]) 
= — 2.4 and —2.5, respectively, and sh ow internal metal- 
hcity spreads of more than 0.5 dex (jKirbv et al.l 120081 : 
lAden et all I2OO90 . Both galaxies contain a popula- 
tion of RR Lyrae stars, implying t he presence of stars 
that are at least as old as 10 Gyr (jMoretti et al.l [20091 : 
iMusella et al.l [2012). Ground-b a,sed imaging further 
suggested old ste llar popula tions (iColeman et al.l l2007t 
Sand et al.l 120101: Td c Jong eFaD 12010: Okamoto ct atl 

llDe 



20121) . IColeman et al- (.2007f ) and iDeason et al- (,201 



suggest that Hercules is tidally disrupting due to its elon- 
gated shape and velocity gradient at large radius. Since 
tidal processes are independent of stellar mass this should 
not affect the IMF analysis, and our HST observations 
described below are well within the gravitationally bound 
region of the object. B ased on the s ame H ST obser- 
vations described below, IBrown et al.l (|2012f ) confirmed 
that stars in these two UFDs are exclusively 11 Gyr or 
older. 

We obtained deep optical images for Hercules and Leo 
IV using the F606W and F814W filters on the HST/ ACS 
Wide Field Camera between August 2011 and January 
2012 (GO-12549, PI: Brown). The total exposure times 
were 25625 and 41060 seconds, respectively. We ob- 
tained two ACS tiles for Hercules and one ACS tile 
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Fig. 1.— The HST/ACS CMDs of two UFD galaxies, Hercules (left) and Leo IV (right). For the IMF analysis, we include stars below 
the sub-giant branch and above the 66% and 75% completeness limits for Hercules and Leo IV, respectively (green points). The axes are 
labeled in both STMAG and Vega magnitudes. The blue line is a representative isochrone of 13.6 Gyr and the mean metallicity of each 
galaxy (Table 1). Blue crosses indicate stellar mass in units of Mq on the main sequence. 



for Leo IV. All images were dithered to mitigate detec- 
tor artifacts and enable resampling of the point spread 
function (PSF). Our image processing includes the lat- 
est pixel-based correction for charge-transfer inefficiency 
(| Anderson &: Bediiil 120101 ). We co-added the images for 
each filter in a given tile usin g the IRAF DRIZZLE pack- 
age (jFruchter fc HooMl2002( ). with masks for cosmic rays 
and hot pixels derived from each image stack, resulting in 
geometrically-correct images with a plate scale of 0.03" 
pixel"^ and an area of approximately 210" x 220". 

We performed both aperture and PS F-fitting phot ome- 
try using the DAOPHOT-II package (|Stetsonlll987h . as- 
suming a spatially-variable PSF constructed from iso- 
lated stars. The final catalog combined aperture pho- 
tometry for bright stars with photometric errors < 0.02 
mag and PSF-fitting photometry for the rest, all nor- 
malized to an infinite aperture. Due to the scarcity 
of bright stars, the uncertainty in the normalization is 
0.02 mag. Our photometry is in the STMAG system: 
m\ = — 2.5x logio/A — 21.1, except where we explicitly 
state Vega magnitudes (Figure [J) . The catalogs were 
cleaned of stars with poor photometry and background 
galaxies. Sources were rejected based on photometric 
errors (< 0.1 mag) and the DAOPHOT x and sharp pa- 
rameters. We also reject stars with bright neighbors and 
those falling within the profiles of extended background 
galaxies because these are noisier than isolated stars of 
comparable magnitude. We apply the same rejection cri- 
teria in determining the completeness of our data de- 
scr ibed below. For mo re details on the data reduction 
see lBrowneFail (|2009( ). 

We performed extensive artificial star tests to evalu- 
ate the photometric scatter and completeness for each 



galaxy. These tests employed the same PSF model and 
PSF-fitting routines used in the construction of the ob- 
served catalogs, including the algorithms for culling poor- 
quality stars, image artifacts and corrections for charge 
transfer efficiency. Artificial stars were inserted into 
the image with appropriate reductions in signal (and 
noise) due to charge transfer inefficiency. Stars were 
then blindly recovered an d corrected for charge transfer 
(| Anderson &: Bedinll2010f) . The artificial stars were in- 
serted over a wide range of color (—1.3 < TTieoe — "m-siA < 
0.5) and magnitude (32 > rrisiA > 16) with most of the 
stars falling near the observed stellar locus and biased to- 
ward fainter magnitudes, thus providing the most fidelity 
in the analysis region. A total of 5,000,000 artificial stars 
were inserted into each image, spread over thousands of 
passes in order to avoid significantly altering the level of 
crowding and the associated photometric scatter. The 
number of stars recovered from all passes sets the com- 
pleteness fractions given in Table 2. For this analysis, 
we include stars fainter than the red giant branch and 
brighter than the magnitude where photometric errors 
approach the main-sequence width and photometric ar- 
tifacts begin to contribute significantly to the catalog 
(green region in Figure [TJ see Table 2). The faint cutoff 
corresponds to the 66% and 75% completeness limits for 
Hercules and Leo IV, respectively. This magnitude re- 
gion corresponds to a stellar mass range of 0.52— 0.77 Mq 
and includes 2380 stars in Hercules and 1054 stars in Leo 
IV. 

2.1. Constructing the Model IMFs 

We analyze our observed luminosity functions by for- 
ward modeling stellar evolutionary tracks with an ana- 
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Fig. 2. — One-sigma confidence contours for the Hercules (green) and Leo IV (blue) UFD galaxies. For the mass range probed by our data 
(0.52 — 0.77 Mq), we plot our results for a log-normal IMF (left) and power-law IMF (right). The best fitting combination of characteristic 
mass (mc) or slope (a) and stellar binary fraction (/binary) a^re shown as plus symbols. For reference, we indicate values of rric and a for 
IMFs in the literature. 



lytic form of the IMF and our observational errors. The 
free parameters in our models are the IMF shape and 
the binary fraction. We also allow the distance and red- 
dening of each UFD galaxy to float within the ori e -signi a 
errors of the values determined in 'Bro wn et al.l ()2012D . 
Because our IMF fits exclude stars brighter than the sub- 
giant branch, the resulting va lues of (m — M)v an d E(B- 
V) are slightly different than iBrown et all ()2012| 1 at the 
hundredth of a magnitude level, providing a better fit to 
the main-sequence region. These values are listed in Ta- 
bl e 1. We fix the gala ctic age to 13.6 Gyr, as determined 
in 'Br own et al.l (|2012[ ) . We further fix the distribution of 
metallicities in each UFD galaxy based on the spectro- 
scopic meta l licity distri bution funct i ons d etermined by 
iKirbv et al.l ([20081 ) and iKirbv et~an (J20ll . We explore 
these assumptions further below. 

We simulate the observations beginning with new 
Victoria- Regina evolutionary tracks (VandenBerg et al., 
in prep.) which assume a solar heavy-element mixture by 
lAsplund et al.l (|2009l ). with 0.4 dex enhancements in the 
a-element abundances, and then scaled to the [Fe/H] val- 
ues of intere st. These are com puted using the same code 
described in iVandenBerg et^ al. (2012). These isochrones 
are transformed into the observed STMAG magnitude 
system using the MARCS model atmosphere library and 
the throughput curves for the ACS F606W and F814W 
filters. The transformation is calibrated to observations 
of globular clusters in the same filters, with 1% agree- 
ment over the main sequence, subgiant branch, and RGB. 

We assume two different functional forms of the IMF: a 
single power-law with free parameter a (where a Salpeter 
IMF is a = 2.35, see Introduction for the functional 
form), and a log- normal model with parameters a and 
characteristic mass ttt-c (where a Chabrier IMF for single 
stars is (T = 0.69 and rUc = 0.2 Mq). In the case of the 
log-normal IMF, our stellar mass range is too small to 
determine the distribution width and, for the purpose of 
this paper, we keep this value fixed at the Chabrier value 



of cr = 0.69. We note that for the Milky Way IMF, cur- 
rent data cannot sig nificantly differentia te between these 
two function forms (jKroupa et al.l[20TTI ). 

The shape of the IMF is somewha t degenerate with th e 
presence of unresolved binary stars (jKroupa et al.|[T99l[ ) . 
We have therefore opted to fit both the IMF shape and 
the binary fraction simultaneously. We define the binary 
fraction, /binary, as the fraction of unresolved systems 
which are binary (e.g., if 25 out of 100 systems are binary, 
then /binary = 0.25 and the total number of individual 
stars is 125). We draw binary companions from the full 
distribution of stellar masses, meaning that the mass ra- 
tio distribution is flat. Because we are explicitly fitting 
for binary fractions, the reported IMF parameters below 
pertain to the single-star IMF, rather than the system 
IMF. 

We convolve the IMF with our known observational 
errors and completeness estimates, resulting in a stel- 
lar probability distribution for a single age and metal- 
licity. While the UFDs likely had a period of star for- 
mation which lasted at least 100 Myr a nd at most 2 Gyr 
(jVargas et al.ll2013t IBrown et al.ll2012| ). our assumption 
of a single age has insignificant effects on the IMF, be- 
cause nearly all of the weight in our fit comes from stars 
below the main-sequence turnoff of the oldest population 
which change by less than 0.01 mag over this time period. 
For the fixed single age of each UFD, we linearly combine 
tracks for different metallicities to recreate the spectro- 
scopic metallic ity distr ibution function (MDF) as deter- 
mined bv Kirbv et all ([2 11) b ased on Keck/DEIMOS 
data from lSimon fc Gehal ()2007[) . The UFDs have signif- 
icant internal metallicity dispersion, over 0.6 and 0.7 dex 
for Hercules and Leo IV, respectively. Since the MDFs 
used in this analysis are based on fewer than 25 stars in 
each galaxy, we roughly estimate the IMF error which 
may be introduced by metallicity errors. We simulate 
two CMDs based on our observational errors and com- 
pleteness values, assuming a Salpeter IMF. For one sim- 



ulated dataset, we use the observed Hercules MDF. For 
the second dataset, we assume a single metallicity of 
[Fe/H]= —2.4. We then recover the IMFs assuming the 
same single metallicity isochrone. The IMF slopes for 
these two simulations are less than 0.1 different in the 
power law slope which is far less than the 1-cr error on 
this parameter determined below. 

The final comparison between model and observations 
is done in CMD space by binning each into Hess dia- 
grams. The Hess diagrams have axes of magnitude and 
color with 0.1 mag bins in each axis. The model and ob- 
served Hess diagrams arc the n compar e d usin g the Max- 
imum Likelihood statistic of iDolpliinl (|2002[ ). using the 
stars defined as members of each galaxy (Figure [T]) . We 
search a grid over two free parameters: the IMF slope 
(a for a power-law IMF) or characteristic mass (rric for 
a log- normal IMF) and the binary fraction /binary The 
best-fit model is that with the minimum Maximum Like- 
lihood statistic over the two free parameters (Figure [3]). 

We do not correct our inferred IMFs for dynamical 
evolution. Dynamical two-body relaxation processes, 
such as mass segregation and evaporation, are mass 
dependent and will alter the initial stellar mass func- 
tion over a relaxation times cale, which we calcu l ate ac - 
cording to Equation 1.38 in iBinnev fc Tremaind ()2008[ ): 
irciax ^ (N/logioN) * {R/(j)- We assume the values for 
the radius (i?) and stellar velocity dispersion {a) using 
the values listed in Table 1. The resulting timescales are 
on the order of 10^^ years, a few hundred times the age of 
the Universe. Thus, two-body relaxation processes can 
safely be ignored in the UFD galaxies. 

3. RESULTS 

The HST/ ACS imaging for Hercules and Leo IV ex- 
tends roughly three magnitudes below the main-sequence 
turnoff (Figure [T|) . Stars included in the IMF analysis 
correspond to stellar masses between 0.52 — 0.77 Af©. 
While this is not a particularly large range in stellar 
mass, it is sufficient to distinguish between various IMFs 
suggested in the literature. 

3.1. Constraints on the IMF 

In Figure [U we plot the one-sigma error contours for 
the binary fraction and IMF parameter a or mc^ in the 
case of the power-law and log-normal IMF, respectively. 
The constraints for Hercules (green) are tighter as com- 
pared to Leo IV (blue) since Hercules is intrinsically more 
luminous and the ACS photometric catalog is better pop- 
ulated along the main sequence. The best-fit solutions 
are consistent between the two galaxies. Based on the 
maximum likelihood statistic, the power law and log- 
normal models provide an equally good fit to the data: 
the log-normal form is slightly preferred but at low (O.Scr) 
significance. 

For a single power law model, we find a best fit slope 
of a = l-2to'.5 for Hercules and a — 1.3 ± 0.8 for Leo 
IV (Figure [21 right panel) . In Figure [3l we plot the ob- 
served luminosity functions for each UFD and overplot 
our best-fitting model (green) compared to a Salpcter 
IMF (a = 2.35, blue) and an extremely bottom-light 
IMF (a = 0.5, red). For Hercules, the IMF slope is shal- 
lower than the Salpeter (a = 2.35) and Kroupa (a — 2.3 
above 0.5 Mq) IMF at the 5.8-ct and 5.4-ct levels. For 



Leo IV, the Salpeter and Kroupa values are 1.9-(T and 
1.7cr from our best fit value. We note that our error 
contours are non-Gaussian. While we strongly rule out 
the Salpeter and Kroupa values for Hercules, we note 
that just below our mass range of 0.52 — 0.77 M0 the 
Milky Way IMF slope turns over. Below 0.5 Mq, the 
Kroupa IMF slope changes to a = 1.3, which is well 
within the UFD values. The exact location of this turn 
over in the Milky Wa.y is un clear. For example the data 
from lBochanski et al.l (|2010l ) plotted in Figure |4] suggests 
a Salpeter-like slope down to 0.3 Mq . We will discuss this 
further in §[5l Our power-law slopes are also shallower 
than recent IMF results in the Small Magellanic Cloud. 
iKahrai et all ()2013f ) found a = I-OOIq;!" over the mass 
range m = 0.37 — 0.93 Mq, which is 2.3-ct away from 
the Hercules IMF. These authors also fit a power law 
only to stars more massive than 0.6 Mq, finding a slope 
a = 2.1 ±0.3. 

We alternatively fit a log-normal function to our UFD 
data, showing the one-sigma confidence intervals in the 
left panel of Figure [2j The best-fitting log-normal 



IMF values are: 

1+2.1 



= 0.4 



+0.9 
-0.3- 



Mq for Hercules and 



mc = 0.415;^ Mq for Leo IV The Milky Way value of 
rUc = 0.2 Mq (|Chabrierll2005| ) is within our 1-cr errors. 
However, the fact that the data prefer a more massive 
characteristic mass is consistent with the results for our 
power law fits, since, in our mass range, a more mas- 
sive characteristic mass provides an overall flatter IMF 
shape as determined for the power law slope. In addi- 
tion, the data for Hercules rule out an extremely 'bottom- 
light' IMF, with nic ^ 2 Mq, as favored bv van Dokkuml 
(|2008D for massive elliptical gala xies, although this result 
has b een subsequently disputed ()van Dokkum fc Conrovl 

[mil). 

In Figure |4j we plot the mass function of our UFDs. 
We compute the mass function by converting the com- 
pleteness corrected F814W luminosity function (Table 2) 
into stellar mass assuming the best fitting single Victoria- 
Regina isochrone and rebinning using constant logarith- 
mic mass bins. While this method does not fully take into 
account metallicity spread or photometric errors (which 
are accounted for in our formal analysis, § 12. ip . this effect 
is small and allows us to compare more directly to liter- 
ature data. Published luminosity functions are available 
for two other galaxies with low mass IMF analyses. We 
convert the luminosity functions for the SMC (Table 1 in 
Kalirai et al. 2012) and the dwarf spheroidal galaxy Ursa 
Minor (Table 9 in Wyse et al. 2002) into mass functions 
in the same manner as the UFDs based on the best-fitting 
single isochrone. Literature data for the mass function 
itself is available only for the Milky Way an d we plot 
data from Table 11 in i Bochanski et all (|2010[) based on 
SDSS photometry of Milky Way field stars. We have ar- 
bitrarily renormalized the mass functions and overplot 
representative analytic fits. While Figure |4| emphasizes 
the limited stellar range over which the UFD IMFs have 
been measured, it visually confirms that the UFD IMF 
slopes are shallower than the observed Milky Way and 
SMC IMFs. We explore the implications of a shallow 
IMF for the UFDs in §|4| and discuss these results in the 
context of other galaxies in §[5| 

3.2. Constraints on the Binary Fractions 
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Fig. 3. — The observed luminosity function for Hercules (left) and Leo IV (right). Errors bars are computed from the observed number 
of stars in each luminosity bin. For comparison, we plot three theoretical power law IMFs, convolved with our observational errors and 
photometric completeness. The fits were normalized to reproduce the number of stars in the observed luminosity function, but here they 
have been normalized at the bright end for clarity. We compare our best-fitting model (green) to a Salpeter IMF (a = 2.35, blue) and an 
extremely bottom-light IMF {a = 0.5, red). 
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Fig. 4. — Stellar mass functions for the five galaxies in which 
the IMF has been measured via direct star counts: the Milky 
Way (blue, Bochanski et al.2010), the SMC (light blue, Kalirai 
et al.2012), Ursa Minor dSph (green; Wyse et al.2002), Leo IV 
(orange; this work) and Hercules (red; this work). Except for Her- 
cules, the vertical normalization is arbitrary. For reference, the 
published power law slopes are shown for each dataset, normalized 
at 0.75 Mq. We note that a power law slope of a = 1 is a flat 
line in this log-log plot. The UFD galaxies show noticeably flatter 
mass functions in this mass range. 



The presence of unresolved binary sta rs can mimic 
a fla ttening o f the IM F at low masses (jKroupa et al.l 
119911; .Bochanski et al.l [201(1) . On the main sequence, 



iKroupa et al.l (|1991t l first demonstrated that unresolved 
binary systems widen the main sequence beyond observa- 
tional errors, with binary systems brighter and redward 
of the single star main sequence. Our best-fitting single- 
power law IMF is degenerate between steeper IMF slopes 
with high binary fractions and shallower (more bottom- 
light) IMF slopes with lower binary fractions (Figure [2]). 
For Hercules, the binary fraction is constrained to be 
/binary = 0.48l^i^ for a single power law or 0.47iJ^:i^ 
for a log-normal IMF. The binary fractions for Leo IV 
are similar, but far less constrained due to the smaller 
number of observed stars (Table 1). While binary stars 
were known to exist in the UFDs based on repeated 
spectroscopic meas urement s of red gi ant branch stars 
(jSimon et al.ll201 it TKoposov et al.ll2011i) . these observa- 
tions did not strongly constrain the binary fraction itself 
(jMartinez et al.l[2011[ ). Our results are the first quanti- 
tative constraints on the binary fraction in the UFDs. 

The binary fractions inferred for the UFDs pertain to 
stars in the mass range 0.52 — 0.77 Mq, predominantly 
K-dwarf stars. In the solar neighborhood, the multi- 
plicity of stars decreases w ith decreasing stellar mass 
(jKraus fc Hillenbrandl I2012J ). with more massive OB- 
typ e stars havi ng binary fractions upwards of /binary = 
0.7 (iPeter et al. 2012), d own to M-stars wi th binary frac- 
tions between 0.2-0.3 (Jans on et al.n 20121. This trend 
may be the result of intrinsic differences in the binary 
fraction as a function of stellar mass, the consequence of a 
single overall binary fraction and random pairings across 
the mass spectrum, or possibly due to binary disruption 
over time combin ed with ages d ifferences between spec- 
tral p opulations (IKroupa et al.1 [1993; Mark s &: K roup^ 
120111) . For K-dwarfs. iDuauennov fc Mavoii (|1991f) sug- 
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Fig. 5. — The power law slope, a, plotted against galaxy velocity dispersion {left), metallicity [Fe/H] {middle) and the observed stellar mass 
range (r i gfef). Data are taken from the same sources as Figure|4] The dotted line in the middle panel is an empirical r elationship suggested by 
IKroupal 1120011) with a zero-point shift to fit these data. Indirect measurements are shown for elliptical galaxies from lConrov iSTvan DokkumI 
1)20121 ). A trend is seen in the sense of shallower, more bottom-light IMF slopes towards less massive and more metal-poor galaxies. 

tion models (e.K.jGovernato et al.ll2012l : lWvithe fc Loebl 
120131: iTevssier et al.ll2013() . particularly for dwarf galax- 
ies. The magnitude of this effect depends in large part 
on the behavior of the IMF outside our observation win- 



gest that roughly half of solar neighb orhood K-dwarf 
systems are binary; a recent analysis by iRaghavan et al.l 
([2010) find a solar neighborhood binary fraction for FGK 
stars of 0.46±0.02. This is in remarkable agreement with 
our UFD binary fractions, despite the significantly lower 
metallicity environment of the UFDs as compared to the 
solar neighborhood. While noting that binary fractions 
as high as /binary = 0.68 are allowed within our 1— cr lim- 
its for Hercules, the best-fitting binary fractions appear 
to be a more 'universal' prop erty than the slope of the 
IMF (IMarks fc Kroupal[2QTl . 

4. IMPLICATIONS OF A SHALLOW IMF IN THE UFDS 

The two UFDs presented have shallower IMF slopes in 
the stellar mass range 0.52 — 0.77 M© than expected from 
a Salpeter or Kroupa IMF and are marginally shallower 
than a Chabrier IMF over the observed mass range. This 
difference is seen visually in Figure [31 For an IMF slope 
a = 2.3, we would have expected to observe 2900 stars as 
compared to the 2380 stars observed for Hercules with an 
inferred IMF slope of a — 1.2. Many properties of the 
UFDs have been calculated assuming a standard IMF. 
Our results have implications for the stellar mass, mass- 
to-light ratio supernova rates and formation models of 
t hese low luminosity systems . 

iMartin et al.l ()2008D calculated the stellar mass of Her- 
cules and Leo IV assuming both a Kroupa and Salpeter 
IMF. Since the Kroupa IMF has a shallower IMF slope 
below 0.5 Af©, the resulting total stellar masses are 50% 
that of a Salpeter IMF. If we instead use our IMF 
to determine stellar mass, assuming that a power law 
slope of a = 1.3 applies over the full mass spectrum, 
our calculated stellar masses have 40% the mass of a 
Salpeter IMF. Alternatively, the best-fit log-normal fit 
with TOc = 0.4 Mq has 44% the stellar mass as a Chabrier 
IMF with rUc = 0.2 Mq. 

Shallower IMF slopes have implications for models of 
UFD galaxy formation. An IMF is assumed in calcu- 
lating the amount of available supernova energy which 
is an important physical process in many galaxy forma- 



dow. If the shallow slope measured at low masses applies 
across the whole mass spectrum, the number of super- 
novae expected per luminous star would increase over 
a Salpeter or Kroupa IMF. Assuming a Milky Way-like 
IMF for the UFDs under-estimates the effects of super- 
nova feedback. 

5. COMPARISONS TO OTHER IMF STUDIES 

The IMF has been measured via direct star counts 
in the stellar mass range 0.5 to 0.8 M© for five distinct 
galaxies: the Milky Way, SMC, Ursa Minor, and the two 
UFDs presented here. As discussed in §[331 we directly 
compare the observed mass functions for these systems 
in Figure [4l We next compare the published analytic fits 
to the IMF, focusing for simplicity only on the power law 
slope (a). 

In Figure [5] we plot the power-law IMF slope as a 
function of the galactic velocity dispersion (a), average 
metallicity ([Fe/H]) and stellar mass range over which 
each measurement was made. The velocity dispersion 
and metallicity of our UFDs are g iven in Table 1. Th e 
SMC IMF results are taken from iKalirai et al.l ()2013D . 
with the galactic v elocity dispersion and metallicity from 
[Harris fc Zaritskvl (2006). I MF results^ for the Ursa Mi- 
nor dwarf galaxy are from lWvse et al . (2002), with the 
galactic velocity dispersion f r om [ Wolf et al. (20101 and 
metal licity from [Kirbv et al] ([20111 ). While IWvse et al.l 
([2002|) do not provide an error bar on the IMF slope, we 
calculate this value by fitting a linear function to the data 
shown in Figure [H For the Milky Way, we use the fitted 
IMF values from Bochanski ct al. ([20101 ) based on a sam- 
ple of Milky Way field dwarfs spanning [Fe/H] = —0.1 
to —0.6. We plot the Milky Way at a velocity dispersion 
of 220/a/2 km s ~\ assuming the disk is embedded in an 
isothermal halo ([Burstein et al.[[l997[ ). 

Figure[5]shows a clear trend in the IMF power law slope 



as a function of galactic velocity dispersion and metal- 
licity. The power law slope becomes increasingly shal- 
low (bottom-light) and less Salpeter-like with decreasing 
galaxy mass/metallicity. This confirms the visual im- 
pression observed in Figure ID As noted in § 13.11 Her- 
cules is 5.4-(T different than the Kroupa IMF and 2.3-a 
away from the SMC value. This is the first clear evidence 
for IMF variations with galactic environments based on 
direct star counts. 

We have not included globular clusters in our com- 
parison in Figure [5l The present-day mass functions of 
many globular clusters have been measured down to 0.3- 
OAMq and these systems are of com p arable metallic- 
ity to the UFD galaxies. iPaust et"an (|2010[ ) and oth- 
ers have measured the present day IMF of Milky Way 
globular clusters, suggesting that globular clusters begin 
with a Milky Way-like IMF and that the observed vari- 
ations are related to dynamical evolution. However, this 
does not fully explain the observed correlations and ad- 
ditional mechanis ms, such as gas-ex pulsion from initially 
compact clusters ([Marks et al.ll2012l ) or orbital properties 
(Kiipper et al. 2013), are needed. Because stars clusters 
do not form within their own dark matter halo, it is 
possible that star formation of these systems is more in- 
fiuenced by their parent galaxy. For these reasons, we 
do not include star clusters in our comparisons with the 
UFDs. 

The trends seen in Figure [S] are qualitatively consis- 
tent with results in elliptical galaxies inferred via in- 
direct methods (i.e., based on integrated galaxy light). 
Correlation of IMF slope with galactic velocity dis- 
persion h as been shown via s p ectroscopic line widths 
_^an Dokk um fc Conrovl 201 it iConrov fc van DokkumI 
[2012; Ferrcras ct al. '2013"), fundamental plane relations 
|Cappcllari ct al. 2012a; Button e t al. 2012) and gravi- 
tational lensing ([Treu et al.l 120101 ). These studies sug- 
gest that the most massive ellipticals have bottom-heavy 
IMFs, with IMF slopes steeper than Salpeter. For ex- 
ample, iCappellari et al.l ([2012b[ ) find a systematic trend 
for a sample of 260 ellipiticals with velocity dispersions 
between 65 - 250km s~^ such that the inferred IMF is 
closer to Kroupa/Chabrier for low velocity dispersion el- 
lipticals, and Salpeter or steeper at highest dispersions. 
In Figure [5] we include for comparison the inferred IMF 
slopes from Conroy & van Dokkum (2012, priv. commu- 
nication) based on 35 elliptical galaxies. These authors 
model the IMF as a broken power-law. We plot their 
'q;2', the IMF slope between 0.5 - 1.0 Mq, a comparable 
range as our UFD observations. The plotted error bars 
are the sample standard deviation. 

A critical question raised by Figure [5] is which physical 
property is responsible for the observed trends with IMF 
slope. While we see a clear trend of increasing (steeper) 
IMF slopes with both increasing galactic velocity disper- 
sion and increasing m etallicity, these two ph ysical prop- 
erties are correlated ([Tremonti et al.l [2004[ ). Metallic- 
ity has often been cited as an expected driver of IMF 
variation, such that lower metallicity gas clouds frag- 
ment into fewer lo w mass protos tars, and thus have a 
more shallow IMF. IKroupal ([20011 ) first noted a possible 
correlation between the IMF slope and metallicity for 
Milky Way clusters and suggested an empirical relation 
of a = 2.3 -I- 0.5 [Fe/H]. The slope of this relationship 
fits the trend seen in the middle panel of Figure [5| and 



we have adjusted the zero-point by 0.2 dex to better fit 
the data. For massive ellipticals, the IMF slope appears 
to correlate more strongly with the alpha-element abun- 
dance [Mg /Fe1 than either [Fe/H] or g alactic velocity 
dispersion ([Conrov fc van DokkumI [20121 ). Since alpha- 
abundances are related to supernova Type II enrichment, 
it is plausible that star formation rates or specific star 
formation density is the c ontrolling physical parameter 
([Weidner fc Kroupal[2005[) . For the UFD galaxies which 
have ~ 10^ Lq a nd formed stars over a period be tween 
100 Myr to 2 Gyr ([Brown et al.|[20TllVargas et al.|[20131 . 
the star formation rates (SFR) are between 10~^ < 
SFR < 10-4 Mq year- ^ The Milky Way has a SF R of 
roughly 1 Moyear"! ([Robitaille fc WhitnevI IMol) and 
the most ma ssive galaxies with t he st eepest inferred 
IMFs by Con rov fc van DokkumI ([20121 ) have inferred 
SFRs up to 100 Moyear"^. Determining the physical 
properties which control the IMF in different galaxies is 
a key question for future studies. 

6. CONCLUSIONS 

We have directly measured the IMF in two UFDs 
via HST/ ACS star counts. Since the UFDs are an 
ancient, metal-poo r and nearly single age population 
([Brown et al.1120121 ) with long two-body relaxation times, 
these objects provide a unique environment in which to 
determine the low-mass stellar IMF. We find a power law 
IMF slope over the stellar mass range 0.52 — 0.77 Mq of 
a = 1.2to:5 for Herc ules and a = 1 .3 ± 0.8 for Leo IV, 
where a = 2.3 for a IKroupal ([2002D IMF in this mass 
regime for the Milky Way. Over our mass range, the 
UFDs exhibit a shallow IMF deficient in low mass stars 
relative to the Milky Way. 

Comparing to other galaxies in which the IMF has 
been measured via directly counting stars, we see 
a trend with galactic velocity dispersion and metal- 
licity. The power law slope becomes increasingly 
shallow (bottom-light) with decreasing galaxy veloc- 
ity dispersion and/or metallicity. This trend is qual- 
itatively consistent with results in elliptical galax- 
ies inferred via indirect me thods (iTr eu et aLJ 12010 ; 
van Dokkum fc Conrovl |201U [C appcUar i et al.l 12012a ; 
Dutton et al.l 120121: IConrov fc van DokkumI 120121 ). The 
combined data provide clear evidence for IMF variations 
with galactic environment. This has significant implica- 
tions for galaxy formation models which often assume 
a Milky Way-like IMF. A galaxy-dependent IMF affects 
estimates of fundamental properties such as supernova 
feedb ack rates, stellar niasses and chemical abundances 
(e.g., iFerre-Mateu et aLll2013f) . A critical question for 
future studies is understanding the physical properties 
which control the IMF as a function of galactic environ- 
ment. 

The UFD IMFs presented here cover a relatively lim- 
ited range in stellar mass. The Milky Way exhibits a 
transition to a shallower IMF slope just below our mass 
limits. Testing whether a similar transition exists in the 
UFDs at 0.5 Mq or lower requires deeper HST imaging 
which can reasonably probe the IMF down to the hy- 
drogen burning limit in the closest UFDs. Alternatively, 
an IMF slope transition could have existed in the UFDs 
at higher stellar mass than our observed range. Since 
stars more massive than 0.77 Mq have long since evolved 



off the main sequence, testing this is difficult, but sig- 
natures may be present i n the detailed chemical abun- 
danc es of individual stars ([Vargas et al.ll2013l : iTsuiimotol 
1201 It ) . Both methods are needed to fully characterize the 
shape of the IMF and thus probe the physical conditions 
of star formation in these extremely old and metal-poor 
systems. 

The details of star formation and the resulting IMF 
have long been predicted to depend on the physical 
properties of the stellar birth cloud (e.g., Larson 2005:; 
Henncbcllc fc Chabricr l200l Krumholz ct al. 2011; 
HoDkinsll2Q12D . It is remarkable that direct evidence for 



IMF variations with galactic environment is only now 
being uncovered, and underlines the utility of the UFD 
galaxies for extending the baseline of physical conditions 
in which the IMF can be directly constrained. 
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TABLE 2 
Luminosity Function for Hercules and Leo IV 
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Note. — The luminosity function and photo 
values. We list values only for the magnitude r 
IMF analysis. 



